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Abstract 

An acoustic metamaterial made of a two-dimensional (2D) periodic array of multi-resonant 
acoustic scatterers is analysed in this paper. The building blocks consist of a combination of 
elastic beams of Low-Density Polyethylene Foam (LDPF) with cavities of known volume. Sound 
inside the structure can excite elastic resonances of the material as well as acoustic resonances in the 
cavities producing several attenuation peaks in the low-frequency range. Due to this behaviour at 
the subwavelength regime we can define this periodic array as a Locally Multi-Resonant Acoustic 
Metamaterial (LMRAM) presenting strong dispersive characteristics of the effective properties 
with subwavelength multi-resonant structural units. The results shown in this paper could be use 
to design effective Sonic Crystal Acoustic Barriers with wide tunable attenuation bands in the 
low-frequency range. 

PACS numbers: 43.20.+g, 43.50.-Gf, 63.20.D- 
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I. INTRODUCTION 



Artificially designed subwavelength electromagnetic materials, denoted metamaterials^-, 
have motivated in last years a great effort to develop both theoretically and experimentally 
the acoustic analogue metamaterial^. Recent works have shown that periodic distributions 
of sound scatterers in a fluid, known as Sonic Crystals (SC)£, can be employed to design 
these acoustic metamaterials. Theoretical predictions and experimental data have shown 
in SC the appearance of frequency ranges, known as Band Gaps (BG), in which sound 
propagation is not allowed^ - -. However, in the regime of large wavelengths, in comparison 
with the separation between the scatterers, SC behave as effective homogeneous acoustic 
medium^. 

The pioneering work of Liu et al.- proposed a novel three-dimensional (3D) periodic 
arrays of coated spheres that exhibited attenuation bands, whose respective wavelength was 
about two orders of magnitude larger than the periodicity of the structure. The origin of 
this phenomenon has been explained as due to the localized resonances associated with each 
scattering unit. Then, sound attenuation within the stop bands increases with the number 
and density of the local resonators, whereas the resonance frequency can be tuned by varying 
their size and geometry. These results introduce the way towards the acoustic analogous of 
the electromagnetic metamaterial, and since the acoustic wavelength of the attenuation band 
in the subwavelenght regime is much longer than the lattice constant of periodic system, one 
can define a Locally Resonant Acoustic Metamaterial (LRAM) whose effective properties 
can provide an accurate and simple description of the wave interaction with the associated 
LRAM. Since the sound speed in the LRAM is proportional to a/ K eff/Peff, where n e ff 
and p e ff are effective modulus and mass density, respectively, a negative K e ///p e // implies 
an exponential wave attenuation (evanescent behavior). Although for natural materials k 
and p are positive numbers to maintain structural stability, for the effective medium with 
low-frequency resonances can appear some differences. It has been shown previously that 
low-frequency band gap was induced by an effective mass density that becomes negative near 
the resonance frequencies, giving rise to exponential attenuation of wave. This is in contrast 
to the acoustic metamaterials formed by SC with no resonant scatterers where the effective 
properties are real and positive. The negative p e ff is seen to result from the coupling of 
travelling waves with the local resonances. In particular, the local resonances enable the 

2 



strong coupling even though the size of the locally resonant unit is much smaller than the 
wavelength in the matrix, in contrast to the Raleigh scattering case. 
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FIG. 1. Experimental data, (a) IL of one U-profile. (b) and (c) show the IL of a SC made of different 
numbers of U-profiles in a triangular array, a =12.7 cm, measured at 0° and 30° respectively. The 
inset shows the transversal view of a U-profile. 



In this work we present a new kind of acoustic metamaterial made of a two-dimensional 
(2D) periodic array of multi-resonant acoustic scatterers. This Locally Mult i- Resonant 
Acoustic Metamaterial (LMRAM) have strong dispersive characteristics of the effective 
properties with subwavelength multi-resonant structural units. The building blocks of this 
acoustic metamaterial consists of a combination of elastic beams of Low-Density Polyethy- 
lene Foam (LDPF) with cavities of known volume. Sound inside the structure can excite 
elastic resonances of the material as well as acoustic resonances in the cavities producing 
several attenuation peaks in the low-frequency range. 

The paper is organized as follows. First, in section [III we characterize the resonances of a 
single 2D elastic beam (EB) as well as the resonances of a rectangular cavity. These results 
allow us to phenomenologically explain the resonances appearing in the scatterers analysed 
in this work. Because of the complex geometry of the scatterer we numerically solve the 
acousto-elastic coupled problems. The numerical model is shown in section IHIl Using this 
model we analyse the behaviour of isolated scatterers as well as the periodic arrays made 
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of them. The experimental validation of all the numerical simulations are shown in section 
IIV1 Due to the resonant behaviour of the periodic structure in the subwavelength regime 
we consider this structure as a Locally Multi-resonant Acoustic Metamaterial in section |Vj 
Finally, the conclusions of the work are presented in section IVII 



A. Motivating results 



The scatteres analysed in this work will be called U-profiles due to the geometrical shape 
of the scatterer (see Figured^). In this Section, the interest is focused on the behaviour of 
the SC made of U-profiles, in the subwavelength regime, meaning that ka << it, where a is 
the periodicity of the array and k is the wavenumber. Hereinafter, this range of frequencies 
will be called low frequency range. 

Figure shows the acoustic response of a commercial scatterer made of the recycled 
material: LDPE closed-cell foam (see inset). To do that we have measured the Insertion 
Loss (IL) defined as the difference between the sound level recorded with and without the 
sample. The red dashed line in the upper graph illustrates the IL of a U-profile. One can 
see the existence of two attenuation peaks appearing in the low frequency range, around 700 
Hz and 1100 Hz. These peaks will be called first and second attenuation peaks, respectively. 
The nature of both can be understood by analysing the eigenvalue and scattering problems 
for the basic geometrical shapes such as rectangular elastic beam and rectangular cavity, as 
it will be seen later. 

On the other hand, this resonant behaviour could be used to improve the acoustic be- 
haviour of the SC introducing attenuation peaks in the range of low frequencies, independent 
of the incidence direction of the wave as in the case seen in the previous Section. The atten- 
uation peaks shown in Figure [T]are obtained with SC made of U-profiles placed in triangular 
array with lattice constant a = 12.7cm. For this lattice constant it is possible to consider 
that the upper bound of the low frequency range corresponds to the first Bragg's frequency 
of that SC with value 1545 Hz. 
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II. PHENOMENOLOGICAL ANALYSIS 



The nature of both attenuation peaks below 1500 Hz can be understood by analysing 
the acoustical properties of basic geometrical shapes like rectangular elastic beams and 
rectangular resonance cavities. In the next two subsections the resonances of both the 
elastic beam (elastic resonances) and the rectangular cavity (cavity resonance) are analysed. 

A. Elastic resonances 

Consider a 2D elastic beam (EB) made of LDPE foam with L length and t width (see the 
schematic view in Figure |2h.); the density of material p. Considering that the EB is fixed 
in one end (black end in Figured), the vibration modes considering temporal harmonic 
dependence can be analysed by means of the following Equation^"— , 

EI 9 ^ = -\ m u 2 Xx) (1) 

where A m = pLt is the linear mass density of the EB, E is Young's modulus, u n is the 
angular frequency of the mode n related to the frequency as co n = 2nv n and I is the second 
moment of inertia. The EI product is known as flexural rigidity. 




length (m> width (m) 

FIG. 2. Eigenfrequencies of an elastic bar of LDPF, density p =100kg/m 5 , Young's modulus 
E =0.35GPa and Poisson's ratio v =0.4. (a) Schematic view of the EB. (b) Dependence of the 
first eigenfrequency on both the length and the width of the EB made of LDPF. 
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The eigenfrequencies of the EB can be obtained from the following Equation, 

cos(fc n L) cosh(£; n L) + 1=0; (2) 



where k n = \Jw\pLtj 'EI. 

From the values of k n and considering that the material of the EB is LDPF, the resonance 
frequencies of the first and second modes are v\ = 693.8 Hz (ki = 28.41 m _1 ) and z/ 2 = 4348 
Hz (&2 = 71.12 m" 1 ) respectively. Using Taylor's series, it is possible to approximate the 
first mode as 



fci 



'12 



(3) 



This first low-frequency solution is particularly interesting for this work. 

In order to validate the use of this analytical approximation, we have compared the results 
obtained for the first resonance peak with those obtained using COMSOL for several values 
of thickness-length ratio {t/L). These values and the relative differences are shown in Table 
[H For the analysed geometry the Euler-Bernoulli equation gives frequency resonance of the 
elastic beam within 1% of the COMSOL results. 

TABLE I. Percentage of the relative differences between analytic and numerical results 
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The resonances of a EB made of a fixed material can be tuned by varying its geometrical 
parameters. Figure [2] shows the dependence of eigenfrequency of the first mode on both the 
length and the width of the EB made of LDPE foam. The black point marks the position of 
the first resonant frequency for the geometrical properties of the EB analysed in this Section 
(t = 0.01 m and L = 0.066 m). One can observe that the bigger the length, the lower the 
frequency of the first mode, and that the bigger the width, the higher the frequency of the 
first mode. The intersection point of the black lines corresponds to the eigenfrequency for 
the EB with the geometrical parameters considered in this work. It can be also observed 
that the eigenfrequencies grow linearly with the width of the EB, however the growth with 
the length is not linear. 



B. Cavity resonances 

Another interesting property of the U-profiles is that they present, in addition to the 
elastic properties of the material, a cavity where sound could be localized due to resonances. 
Several works in the literature have analysed the effect of cavity resonators in periodic 
structures. In addition to the BG of the periodicity, the systems made of resonators show 
low frequency attenuation bands produced by the resonances of Helmholtz or split-ring 
resonators. Due to both the control of the resonances of the EB and the resonance of the 
cavity, periodic structures made of U-profiles elastic scatterers can be easily tunable in the 
range of low frequencies. 



(a) 



B, 




FIG. 3. U-profile elastic scatterers. calculated using FEM. (a) Geometrical shape and parame- 
ters characterizing the size of the U-profile. (b) Acoustic field inside the cavity for the resonant 
frequency v r . 
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The eigenfrequencies of a rectangular cavity with several boundary conditions has been 
widely analysed in the literature. As one can see in Figure [3^, the cavity of the U-profile 
has a length and a width equal to l x = 0.066 m and L 2 = 0.04 m respectively. To solve the 
problem, one can consider that the walls of the U-profile are perfectly rigid. Thus, Neumann 
boundary conditions should be considered at the boundaries, and Dirichlet conditions in 
the boundary in the open side of the U-profile. The solution of the analytical problem 
results in an eigenfrequency problem whose fundamental mode has a eigenfrequency equal 
to v r = c a i r /{Al x ). However, we have shown that the air immediately outside the end of 
the cavity takes part in the acoustic oscillation. This air makes the cavity appear to be 
acoustically somewhat longer than its physical length. This effective length gives rise to a 
shift of the resonance frequency, for this reason it is called in the literature, the end correction 
of the cavity. In order to compute the correct resonance frequency, this effective length and 
the corresponding frequency correction have to be considered. A more rigorous analysis of 
the cavity would be required to find the exact resonance frequencies, but it has been assumed 
that the resonance is affected by the end correction of the cavity of the U-profile following 
the next equation^ 



is the end correction used 

inM. 

The adaptation of the equation C29 in reference^ can be 
seen in the |A] Then, in the case of the cavity considered in this work, the frequency of the 
first mode is v T =1155 Hz. 

III. NUMERICAL RESULTS 
A. FEM model 

We start this Section by analysing the propagation of acoustic waves inside periodic 
structures made of solid scatterers, B, embedded in a fluid host, A, using FEM. Due to 
the physical properties of the host material, the eigenmodes of the whole system are pure 



c, 



(4) 



A{l x + 5Y 



where 
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longitudinal waves, while transverse modes cannot propagate. Then, the governing Equation 
in A is 

where p is the pressure, pa is the density and ca is the sound velocity in the host material. 

The propagation of elastic waves inside the scatterers, locally isotropic medium, is gov- 
erned by 

where ps is the density of the elastic material and ui is the ith component of the displacement 
vector. The stress tensor is defined by 



<y%j = ^BUuSij + 2p B v>ij 

1 J dui duj 

2 \ dxj dxi 

where As and pb are the Lame coefficients. 

In this problem the acoustic wave is incident on the scatterer and then the pressure acts 
as a load on the elastic medium. On the other hand, the elastic waves in the scatterer act 
as an additional acceleration on the acoustic field. In order to simultaneously solve [7] and M 
we introduce the following boundary conditions, 

g^ldB = Paoj un (9) 
(?ijnj\dB = ~V n i- 

where OB is the boundary of the medium B and n is the outward-pointing unit normal 
vector seen from inside the scatterer medium. 

To solve the stated problem COMSOL MULTIPHYSICS has been used, as well as a 
finite-element analysis and solving software package. FEM is a good technique when the 
shape of the involved objects is complicated and several physical problems are coupled. 
In the numerical problem, the domain in which the solution is obtained was surrounded 
by Perfectly Matched Layers (PML) region in order to emulate the Sommerfeld radiation 
condition in the numerical solution 15 . 
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B. Scattering problem 



1. Single scatterer 

First of all, the frequency response of the U-profile will be analysed. The geometry of 
the U-profiles was implemented using the CAD tools of COMSOL. In the inset of Figure 
HI one can see the considered model of the scatterer, this geometry slightly differs from the 
real one. In the numerical model, a plane wave impinging the scatterer from the left side 
has been considered and the IL at a point behind the scatterer is calculated. 

Figure H] shows the numerical results obtained using COMSOL. One can observe similar 
frequency response to the one experimentally observed in Figure [H Two attenuation peaks 
numerically calculated appear near 700 Hz and 1200 Hz. A sensibility analysis varying 
the geometry of the U-profile was done, and one can observed that the first peak is very 
sensitive to changes in geometry. However the second one does not substantially change 
with the variations of the geometry of the U-profile. These changes are in good agreement 
with the predictions of the resonances of a rectangular elastic beam profile with a resonant 
cavity. 
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FIG. 4. Numerical results of single scatterer. IL produced by a U-profile. The inset shows a image 
of the numerically modelled U-profile. 

As we have previously explained, the first peak corresponds to the resonances of the EB 
of the U-profile, therefore small changes in the geometry can produce high changes in the 
resonant frequencies (see Figure [2]). This provides a powerful design tool in the first of the 
low frequency attenuation peaks. 
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On the other hand, if we compare Figures [T] and H] we can observe again a difference 
between the numerical and experimental frequencies of the second peak. The modelling of 
wave scattering in FEM should be able to describe any effects within the theory of linear 
acoustic and solid mechanics. Use of the end correction could be only appropriate for 
analytical approximation which neglect the interactions between the cavity interior and the 
environment. Any discrepancies observed between numerical results and experiments are 
most likely due to the idealization of the scatterer material. 



C. Infinite periodic array of U-profiles 

We analyse the propagating properties of a periodic arrangement of U-profiles by means 
of the application of the Bloch periodic boundary conditions in a unit cell. That means 
studying the band structures for a periodic array of U-profiles. Figure |5] shows the Band 
structures of a periodic arrangement of U-profiles arranged in triangular lattice of a =12.7 
cm. The black line represents the band structures considering the perfectly rigid U-profiles, 
this means that we can consider Neumann boundary conditions in the wall of the scatterer. 
One can observe that this arrangement presents the pseudogaps related to the periodicity 
as well as the stop band due to the resonance of the cavity. 

If the elastic properties of the U-profiles are considered, then the band structures are 
represented by the blue continuous line. One can observe the BG due to the periodicity 
(~1600 Hz), the stop band of the resonance of the cavity (1100 Hz) and the resonance of the 
elastic beams of the U-profile (700 Hz). In the representation the non propagating ranges 
of frequencies are presented by the black surfaces. 

In order to compare the numerical results with those experimentally obtained a new plot 
was added in the right panel, where the IL of a triangular periodic distribution of U-proflies 
was measured in the main symmetry directions: 0° (blue line) and 30° (red line). One 
can observe the low dependence on the direction of incidence of the resonance due to the 
vibration of the elastic beams and the resonance of the cavity, but the directionality appears 
in the frequencies of the BG due to the periodicity. 

11 




X r J X 5 10 15 20 25 30 35 

k Acoustic Attenuation - IL (dB) 



FIG. 5. Band Structures for a periodic arrangement of U-profiles in triangular lattice of a =12.7 
cm. Left panel: Black dashed line represents the band structures for a rigid U-profiles, whereas 
the blue line represents the bands for the elastic U-profile. The black surfaces indicate the non 
propagating ranges of frequencies. The Brillouin zone is defined as V = (0,0), X = (tt, — tt/\/3), 
J = (4-7r/3,0). Right panel: Measured IL of a triangular lattice of elastic U-profiles measured in 
the two main symmetry directions, 0° (blue line) and 30° (red dahsed line) 

1. Finite periodic array of U-profiles 

Once the numerical results of the acoustical behaviour of an isolated U-profile have been 
analysed, the next step is to analyse a periodic distribution of this elastic scatterers following 
a triangular lattice with lattice constant a = 0.127 m. Here, a plane wave impinging from 
the left side is considered and the numerical domain is again surrounded by a PML region. 
Thus the numerical solution accomplishes an approximated Sommerfeld condition. 

The blue continuous line in Figure [6] shows the numerically predicted IL of a finite struc- 
ture made of 6 rows of 10 U-profiles for an incident wave in the direction of 0°. By comparing 
the results of the scattering of an isolated U-profiles (see Figure [T]), one can observe that 
the resonance of the elastic beams, as well as the cavity resonance, have been increased due 
to the increase in the number of resonators. Moreover, an attenuation peak around 1600 Hz 
appears and it can be related to band gap of the array. 

Open circles in Figure [6] shows the measured IL for the same array as the numerically 
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calculated is also plotted. One can observe a good correspondence between the attenuation 
peaks numerically predicted with those experimentally obtained. However, the experimental 
attenuation peak related to the elastic resonances presents a higher attenuation level than 
the one numerically predicted. A possible explanation for this effect could be the existence 
of some absorption effect of the material that it is not considered in the model. 




400 600 800 1000 1200 1400 1600 1800 2000 2200 

Frequency (Hz) 



FIG. 6. IL of an array of U-profiles. Blue continuous line represents the numerically predicted IL 
and open circles represent the measured IL. A plane wave impinges the structure from the left side 
to the right side. 



Then, the resonance effect of the scatterers is not destroyed by the multiple scattering 
inside the structure, therefore one can combine resonances with multiple scattering in order 
to obtain several attenuation peaks. As we explained in the beginning of this Section, 
an interesting feature of these scatterers is that they show two resonances in the range of 
frequencies below the first BG of the periodic structure. 

IV. EXPERIMENTAL RESULTS 

The acoustic attenuation capabilities of the single LDPE foam scatterers as well as of 
the periodic structures of these scatterers have been measured in terms of IL. The sample 
is excited by white noise. 
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A. Single scatterer 

The experimental analysis of the single LDPF scatterers has been divided into two parts, 
acoustic and vibrational analysis. 

6— 

4 

m 

S 2 
-j 

o 

-2 
006 

u> 

£0.04 

c 
o 

jS 0.02 
<D I 

ai U: 
o 'J 

< 0/ 



FIG. 7. Experimental results of a single scatterer. (a) IL (dB) measured behind the scatterer. Red 
dashed line represents the experimental results and Blue line represents the numerical simulation 
using FEM, (b) Experimental measurements of the vibration of the EB of the LDPE foam scatterer. 
Blue line represents 0° of incidence and red dashed line represents 30° of incidence, (c) Sound level 
map measured inside the cavity for the resonant frequency for v =1104 Hz. Step Ax = Ay =1 cm. 

The IL of the single scatterer is shown in Figure Eh in red dashed line. The real shape 
of the scatterer was modelled using FEM. Using the acoustic-elastic coupling previously 
presented, the IL was numerically obtained (see the blue line in Figure Ek). One can observe 
that the numerical method previously used is in good agreement with the experimental 
results. The observed discrepancy between the numerical and experimental results in the 
second peak can be explained by the lack of precision in the representation of the scatterer 
profile in the numerical method. 

In Figure [T^i one can see the presence of the cavity and the EB resonances. In order to 
experimentally explain these resonances, we have measured on the one hand the vibration 
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of the EB with an accelerometer, and on the other hand, the acoustic field inside the cavity 
of the LDPE foam for the resonant frequency of the cavity. In Figure [7b, one can observe 
the vibration of the EB for two different incidence direction of the acoustic wave. Blue 
continuous line represents the vibrations of the wall of the U-profile for the incident wave in 
the TX direction. Red dashed line represents the vibrations of the walls for a wave impinging 
in the TJ direction. The vibration of the wall of the elastic beam increases at the resonant 
frequencies. Figure [7b shows the increasing in the vibration of the wall in the resonance 
of the material and in the cavity resonance independently of the incident direction of the 
acoustic wave. 

Figure shows the acoustic field inside the cavity of the LDPE foam obtained by moving 
the microphone with our robotized acoustic measurement system in 1 cm steps inside the 
cavity. The field inside the cavity is similar to the one numerically obtained in Figure [3H. 
The resonance of the cavity induces the vibration of the walls as observed in Figure [7b. 



B. Periodic array 



1. Dependence on the number of resonators 




FIG. 8. Experimental measurement of the IL for determining the dependence of the attenuation 
peaks on the number of scatterers. Open colored circles represent the IL for six structures made of 
different number of rows (from 1 to 6 rows of 10 cylinders per row). IL measured 1 m away from 
the end of the complete structure, (a) Measurement in the TX direction, (b) Measurements in the 
LJ direction. 



One can expect that both the resonant effect and the multiple scattering depend on the 
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number of scatterers in the array. To prove it, we built and measured six configurations 
with an increasing number of scatterers. The final structure presents 6 rows of 10 scatterers 
per row. We measured the IL at the same point for the six structures that have from 1 to 
6 rows respectively. Figure E] shows these experimental results. The coloured open circles 
represent the IL measured 1 m away from the end of the complete crystal. Blue open circles 
show the IL of a structure made of 1 row of U-profiles for 0° incidence whereas the black 
open circles represent the IL of structure made of 6 rows of U-profiles for 30° of incidence. 

We can observe that both resonances due to the elastic material and the cavity, depend 
on the number of cylinders in the structure. Also, in the case of only 1 row, where there is 
no periodicity in 2D, the resonance peaks are present in the attenuation spectrum whereas 
Bragg's peaks do not appear. The attenuation spectra of structures made of rigid scatterers 
always present ranges of frequencies where there is sound reinforcement, meaning negative 
IL. However, it is interesting to note that these structures do not present ranges of frequencies 
with sound reinforcement. 

2. Dependence on the incidence direction 

One of the main characteristics of the stop bands produced by periodic arrays is their de- 
pendence on the incident direction. In periodic systems the BG results from the intersection 
of the frequency ranges of the pseudogaps in the main symmetry directions, the upper and 
lower bands of each main direction being dependent on the incident direction. However, it 
is known that the resonance effect must be independent from the incidence direction. 

Here, we measure this dependence of both the resonance and the multiple scattering in a 
periodic array of U-profiles in the direction of incidence. We have especially measured the 
IL of a complete structure for several incident directions, between the two main symmetry 
directions (0° and 30°). Figure |9] shows these experimental results measured in the anechoic 
chamber. 

One can observe in Figure [9] the low dependence on incidence direction of the attenuation 
peaks produced by the resonances of the elastic walls and by the cavity of the U-profile. 
However, one can see that the behaviour of the attenuation peak produced by the multiple 
scattering in the periodic system is highly dependent on the incidence direction. 
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FIG. 9. Experimental measurement of the IL to determine the dependence of the attenuation peaks 
on the direction of incidence of the wave. Open coloured circles represent the IL for four different 
directions, 0°, 10°, 20° and 30°. IL measured lm away from the end of the complete structure (6 
rows of 10 cylinders per row). 



V. LOCALLY MULTI-RESONANT ACOUSTIC METAMATERIAL 

In LRAM, the sound speed is proportional to \/^eff/Peff, where n e ff and p e ff are the 
effective modulus and the mass density, respectively. Depending on the values of these 
parameters the metamaterial presents several responses in frequency. SC made of rigid 
scatterers with no resonant properties can be analysed as an acoustic metamaterials showing 
real and positive effective properties^. However, some interesting differences can appear due 
to the effective medium with low-frequency resonances. 

In order to have a propagating plane wave inside the medium, we should have either 
both positive n e ff and p e ff or both negative n e ff and p e //. Moreover, with these values the 
Poynting vector for a propagating plane wave is defined by 

S = ^ p Vp*= 1 -^. (10) 

If n e ff and p e ff are positive, the Poynting vector, S, presents the same direction as k and 
the Snell law is normally accomplished. However if n e ff and p e ff are negative, S and k 
present opposite directions, and the metamaterial behaves as a left handed material, where 
the negative refraction appears. Physically, the negativity of n e ff and p e ff means that the 
medium displays an anomalous response at some frequencies such that it expands upon 
compression (negative bulk modulus) and moves to the left when being pushed to the right 
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(negative density) at the same time. These materials present unique properties due to the 
double-negative medium, such as negative refractive index and subwavelength focusing^. 

However, if only one of both parameters K e ff and p e // are negative, the sound velocity is 
complex, and the vector presents a complex value. Thus, when the real component of the 
expression of the Poynting vector is negative and sufficiently large, we can observe a narrow 
frequency range, corresponding to the region of negative modulus, where Re(kS) < 0. A 
direct consequence of such behaviour is the exponential wave attenuation in such frequencies. 
It has previously been shown that low-frequency attenuation bands can be induced by an 
effective bulk modulus that becomes negative near the resonance frequencies, giving rise to 
exponential decay of modes 5 . 

In the system studied here, the stop bands at low frequencies are independent of the 
angle of incidence and of the lattice constant 7 . Moreover, there is not transmission wave 
in the resonance frequencies, consequently we do not observe any negative refraction or 
subwalength imaging near the resonance. These properties would imply propagation in 
some frequency region. Thus, we can conclude that, as in the case of reference^, the periodic 
structure made of U-profile could be represented by one of the effective parameters being 
possitive and the other one negative in the resonant frequencies. 

A rigorous parameter retrieval procedure on the line of those developed for electromag- 
netic and acoustic cases will be required to be implemented on this system to obtain the K e ff 
and p e ff- But since the LRAM structure, in our case, is not a far subwavelength in size of the 
operating frequency (A ~ a), such a homogenization of all properties via effective medium 
parameters is difficult. However, one can follow the formalism of the electromagnetic mate- 
rial to phenomenologically analyse the behaviour of the system in the subwavelength regime 
analogously as in reference 5 -. 

Up to the best of our knowledge, the negativity of both the bulk modulus and the effective 
density in acoustic metamaterials is related to the monopolar and the dipolar resonances 
of the building blocks of the metamaterial respectively. The monopolar Helmholtz^ reso- 
nances have been used to design acoustic metamaterials with negative bulk modulus and 
the dipolar resonances of spheres have been used to design double-negative metamaterials^, 
i.e., negative bulk modulus and dynamic mass density. In this work we have shown evidence 
about the different origin of the resonances in the U-profile system being the two attenuation 
peaks produced by monopolar resonances. Therefore both resonances could be considered 
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in the expression of the bulk modulus. 

The acoustic properties of a 2D SC can be mapped into an electromagnetic counterpart, 
where p, v, p, k correspond to H z , E, e, /x, respectively. Following the formalism of the 
electromagnetic metamaterials, one can consider that the systems behave as a metamaterial 
with an effective bulk modulus K e ff(u) in the form, 



K eff [U) 



? U(i 



E 



2u) 



-i 



1 - 



u 2 — Uqj + iTu 



(11) 



where F is the filling fraction, u^j represents the resonant frequencies of the LDPF scatterer, 
T is the dissipation loss in the resonating elements and N res is the number of resonances 
of the scatterers. The loss term F cannot be determined a priori therefore it should be 
experimentally determined. In our system N res = 2 and the resonances are represented by: 

VT2 fpLtY 1/2 
w i = ^tt- I — J (12) 

(13) 



L 2 \ EI 

^02 



2^rc a ^ r 



4(/, + 5)' 

In Figure [TUh . one can observe the effective bulk modulus of the material. The imaginary 
part of the effective bulk modulus presents this particular frequency-dependent response 
which is essential to the range of frequencies where a stop band is expected. 




1 20 10 

Inserion Loss (dB) 



(b) 



(d) 



FIG. 10. Effective parameters, (a) Effective bulk modulus. Imaginary part is plotted in red dashed 
line whereas real part is plotted in blue continuous line, (b) Dispersion relation, (c) Transmission 
coefficient for a slab of metamaterial with L e *f =0.549 m. (d) Measured IL of an array of U-profiles. 
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To obtain the transmission coefficient of a slab of the metamaterial with the bulk modulus 
shown in Equation [TTJ it is necessary to determine the size of the effective thickness of the 
material. The filling fraction of the structure is, 

f = hif^L (14) 

A eff 

where N is the number of scatterers, A cy \ is the area of each scatterer and A e ff is the area 
occupied by the homogeneous scatterer. For a homogeneous scatterer of N scatterers of 
area A cy i in a lattice whose unit cell has an area A uc , the filling fraction gives the following 
Equation 

A-cyl _ NA cy i 
A uc A e ff 

For the homogenized system with the rectangular shape considered in this work (6 rows of 
10 scatterers with triangular periodicity), the effective thickness L e ff follows 

_ \/3a 



J eff 



(16) 



2 ' 

where a is the lattice constant of the square array of the inner structure of the homogeneous 
material. For the parameters considered in this work, L e fj =0.549 m. 

Finally, the density of the effective medium is considered here as the dynamic mass 
density, 

Peff = Ph + p s - f(ph - Ps) ^ 
Ph Ph + Ps + f(ph ~ Ps) ' 

where, p s is the density of the scatterer and ph is the density of the host material. 

It would be interesting to know what is the dispersion relation corresponding to this 

medium with negative elastic modulus. In the regime of low frequencies, the real (x = Re{K)) 

and imaginary (y = —Im(K,)) parts of the bulk modulus can be related to the propagation 

constant of the media as£, 



Re(k) = (V^T? - x) V2 , (18) 

Im(k) = (v^T? + x) V2 . (19) 

In Figure [TUb . one can observe the dispersion relation obtained from the effective elastic 
modulus, K e ff in Equation [TTl using Equations [181 and [T9l Two spectral bands of no prop- 
agating modes are expected in the vicinity of the resonances of the local resonators. 

20 



From the Fresnel Equation of stratified media, it is possible to calculate the transmission 
coefficient of a slab of 0.549 m of the acoustic metamaterial analysed here. In Figure [TUb. we 
can observe the absolute value of the transmission coefficient. One can observe a reduction 
of the transmission around the resonant frequencies. 

In Figure [TOH . we show data from measurements of the IL of a periodic array of U- 
profiles distributed in a triangular lattice. We observe that the attenuation peaks predicted 
by using the effective medium approximation are also experimentally seen. On the other 
hand, we would like to note that the diffraction limit is near to the second peak, and in the 
experimental measurements, one can observe Bragg's peak due to the periodicity. 



VI. CONCLUSIONS 



In summary, the resonances of a scatterer with complex geometry have been studied 
from the analysis of the resonances of simpler geometries. The easy way to design the 
resonances of the scatterers presented in this work opens several possibilities to create an 
arrangement of resonant scatterers that attenuate a wide range of frequencies below the BG 
of the structure. The first possibility is to analyse the acoustical behaviour of a periodic 
array made of several scatterers with different length and cavities. Following the rules for the 
design of the resonance frequency of both elastic beam and cavity resonances, it is possible 
to design tunable stop bands in the propagating range of a SC if it were made of rigid 
scatterers. The results shown in this paper could be used to design effective Sonic Crystal 
Acoustic Barriers with wide tunable attenuation band in the low-frequency range. 
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Appendix A: End correction formula 



The end correction obtained by Norris et al. in reference^ for cylindrical rigid resonators 

is: 

r 2aa i f 2 \ 

5~ log - , Al 

7r \a J 

where a is the angle of the aperture of the cylindrical Helmholtz resonator and a the radius of 
the cavity Taking into account that the resonant frequency of the 2D Helmholtz resonators 
depend on the area and the aperture of the cavity, we have adapted this formula to our case 
considering a cylindrical resonator with the same area and aperture as our U-profile. Then 
the equivalent radius and aperture angle follow 



lxU (A2) 



7T 



a e = 2 arcsin I — — ) . (A3) 



2a e 

These definitions need additional conditions that make the aperture angle a e real. The 

condition | sin(o; e /2)| < 1 that gives L 2 < 4:I x /tt. With this, the end correction follows 

2a e a e ( 2 \ , . 

5 ~ log — . (A4) 



7r \ a 
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